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If the Examiner believes a telephone conference would expedite 
prosecution of this application, please telephone the undersigned at 650-326-2400. 
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VERSION WITH MARKINGS TO SHOW CHANGES MADE 

IN THE SPECIFICATION: 

The paragraph beginning at page 10, line 30 has been amended as 

follows: 

An asterisk indicates the transpose: H\k) , for example. A tilde over a vector indicates a 
3 by 3 skew-symmetric cross product matrix: vw=vx w . is an i by i identity matrix, 
and 0,. is a zero vector of length i and 0. is an i by / zero matrix. 



The paragraph beginning at page 23, line 18 has been amended as follows: 

The implicit Euler integration method is illustrated in the flow chart of 
Fig. 6 for the vector function y = f(y,t) (where y = (q,u) y q representing the position 
states and u the velocity states of the molecular system). The function / includes both 
the multibody system dynamics and the forces such as electrostatic attraction and 
repulsion, van der Waal's forces, and solvation forces. After an entry step 79, the first 
operation step 80 updates the Iteration matrix G . For all implicit integration methods, 
the Iteration matrix G has the form G = I-aJ , where / is the identity matrix, a is 
some scalar function of the timestep h n , the timestep between time t n and t n _ x , and J , 

the Jacobian given by 7=— . For the implicit Euler method, a = h n . In passing, for 

dy 

additional savings in computer time, it should be noted that a very efficient method of 
computing Jacobian matrices from the residual form of equations is covered in previously 
cited co-pending U.S. Patent Appln. No. 10/053,348, entitled "METHOD FOR 
ANALYTICTGTJXCOBIANCOMP 

of even date and is assigned to the present assignee. As in the case of the present 
invention, the same referenced patent application also describes the use of internal 
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coordinates to describe the state of the molecular system. For example, the rotation of 
one part of the molecule is described with respect to another part, rather than with respect 
to an external referenced coordinates. This further increases computing efficiency. 
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